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ABSTRACT 


The "PPN formalism" — which encompasses the post-Newtonian 
limit of nearly every metric theory of gravity — is used to 
analyze stellar stability. This analysis enables one to infer, 
for any given gravitation theory, the extent to which post- 
Newtonian effects induce instabilities in white dwarfs, in 
neutron stars, and in supermassive stars. It also reveals 
the extent to which our current empirical knowledge of post- 
Newtonian gravity (based on solar-system experiments) actually 
guarantees that relativistic instabilities exist. In parti- 
cular, it shows that: (i) for "conservative theories of gra- 

vity", current solar-system experiments guarantee that the 
critical adiabatic index, r j , for the stability of stars 
against radial pulsations exceeds the Newtonian value of 4/3: 

r crit = 4/3 + K M/R , K positive and of order unity; 

(ii) for "nonconservative theories", current experiments do 

not permit any firm conclusion about the sign of r „ - 4/3: 

crit ' ’ 

(iii) in the PPN approximation to every metric theory, the 
standard Schwarzschild criterion for convection is valid. 
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I. INTRODUCTION AND SUMMARY 


Relativistic corrections to Newtonian gravity should induce dynamical 
instabilities in stars with adiabatic indices slightly greater than U/3. 

This fact was first discovered, within the framework of General Relativity 
(GR), by Chandrasekhar (l96ka, b) and independently by Feynman [unpublished, 
but quoted in Fowler (196^)]. More recently Nutku (1969) has shown that 
the same type of instability is predicted by the Brans-Dicke theory of 
gravity (BDT), but that it is slightly weaker (stars are slightly more 
stable) than in GR. If the dynamical relativistic instability actually 
exists, as predicted by GR and BDT, then it plays a fundamental role in 
white dwarfs, in neutron stars, and in supermassive stars [see e.g. Thorne 
(1967) or Zel'dovich and Novikov (1971) for a review]. 

But it is conceivable that neither GR nor BDT is the correct relati- 
vistic theory of gravity. If so, might the relativistic instability not 
exist? Is it conceivable that relativistic effects would stabilize stars 
rather than destabilize them? William A. Fowler has asked this question 
of gravitation theorists so often at Caltech, that we have felt compelled 
to seek a firm answer. The most firm answer possible is one which relies 
heavily on experimental tests of relativistic gravitational effects, while 
assuming nothing (or almost nothing) about which relativistic theory of 
gravity is correct. 

Of course, one cannot work in a complete theoretical vacuum. A minimal 
amount of theory is required to link the relativistic instability in stars 
to solar-system measurements of perihelion shift, light deflection, radar 
time delay, etc. That the amount of theory needed is small, however, one 
can see heuristically by noticing that both the perihelion shift and the 
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relativistic instability are caused by a relativistic strengthening of 
Newtonian gravitational forces. [Stronger gravity than predicted by Newton 
when a star contracts means greater force to pull the star on inward, i.e. 
means less stability; stronger gravity than predicted by Newton when a 
planet approaches close to the sun (perihelion) means greater force to 
"whip" the planet around, and a resultant advance of its perihelion.] 

The purpose of this paper is to derive a quantitative measure of the 
extent to which solar-system experiments imply the existence of the dynami- 
cal relativistic instabilities in stars. The "minimal amount of theory" 
to be used in the derivation is the Parametrized Post-Newtonian ("PPN") 
Framework of Nordtvedt and Will (Will and Nordtvedt 1972; Will 1971a; 


Nordtvedt 1968). 

The PPN Framework is a post- Newtonian theory of gravity with adjustable 
parameters. In Will's fluid version, it has nine PPN parameters, 7, 8, 


a , a , £ , £ , £ , and . The parameter 7 measured curvature of the space- 
2 3 12 3 * 

geometry; 8 measures the non-linearity of gravity; Qf p a 2 , and a ^ measure 


"preferred- frame" effects; £ p £ 2 , £ 3 , and ^ measure the effects resulting 
from a breakdown of conservation laws. For theories which have no "preferred- 


frame" effects, all a' s vanish (Nordtvedt and Will 1972; Will 1971b). For 


theories which have conservation laws for energy, momentum, angular momentum. 


and center- of-mass motion ("conservative thoeries"), all a' s and £'s vanish 
(Will 1971b). The post-Newtonian limit of every "metric theory of gravity" 


Metric theories of gravity are a wide class of theories including (i) 
every theory that satisfies the equivalence principle (laws of physics in 
local Lorentz frames the same as in special relativity), and (ii) every 
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theory that the Caltech group has thus far examined and found to be complete, 
self-consistent, and in agreement with experiment. See Thorne, Will, and Ni 
(1971); Ni (1972) and Will (1972b) for full discussions. 


known to us [except Whitehead's theory which is non-viable (Will 1971c] is 
a special case of the PPN Framework, corresponding to particular values of 
the PPN parameters. Ni (1972) has calculated the values of the parameters 
for a variety of theories, including general relativity, the scalar- tensor 
theories of Bergmann-Wagoner, Nordtvedt, and Brans-Dicke- Jordan, the 
conformally- flat theories of Witrow-Morduch, Littlewood-Bergmann and Nordstrom, 
and the stratified theories of Einstein, Witrow-Morduch, Page-Tupper, Yilmaz, 
Papapetrou and Rosen. 

Experiments to date have placed the following limits on the PPN para- 
meters [see Thorne, Will, and Ni (1971) or Will (1972b) for detailed dis- 
cussion; see also Nordtvedt and Will (1972)]: 


7 = 1.0k ± 0.08 [time delay and deflection experiments except 
that of Sramek (1972)] 

7 = 0.88 ± 0.12 [ Sramek' s (1972) deflection experiment] 

8 = 1.14 [perihelion shift plus time delay experiments] 

[£ _ — ~ \ ol z \ ~ 0,11 [ Kreuzer (1966, 1968 ^measurement 

^ of m . . /m , ] 

active' passive 

|t| 3 | < 0.05 [Kreuzer (1966, 1968) measurement of m active / m passiv£ J 

|0L I <0.2 [Earth rotation rate experiments (Nordtvedt and 
1 Will 1972)] 


( 1 ) 

( 2 ) 

(3) 

00 

2 

( 6 ) 


( 5 ) 


3 


a 2 \ < 0.03 [Earth-tide measurements (Will 1971b)] 


(?) 


|a | $ 2 X 10 ^ [perihelion shift observations (Nordtvedt 
5 and Will 1972)] (8) 
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Kreuzer's (1966) analysis of his data was completely correct, despite a 
recent claim to the contrary by Gilvarry and Muller (1972). Gilvarry 
and Muller err in making a quadratic fit to Kreuzer's data, rather than 
restricting themselves to a linear fit as did Kreuzer. Kreuzer measured 
the expansion of his liquid relative to teflon over a wide temperature 
range and thereby showed experimentally that the quadratic correction to 
the linear behavior must be negligibly small over the small temperature 
range used for the experiment. Moreover, the magnitude of the quadratic 
term obtained by Gilvarry and Muller using their least- squares fits is 
ridiculously large for any but pathological materials. We thank R. H. Dicke 
for a helpful discussion of these points. 


In this paper it is shown that for conservative theories of gravity 
current experimental limits on the PFN parameters — based on perihelion 
shift, time delay, and deflection experiments — guarantee the existence of 
the dynamical relativistic instability in stars; while for non-conservative 
theories the present, experimentally undetermined state of the two PPN 
parameters £ and makes it uncertain whether relativistic effects will 
actually stabilize or destabilize stars. In quantitative terms, a non- 
rotating spherically symmetric star with adiabatic index 

r, -(&H) - (t) (9 > 

8 S 

k 




constant throughout its interior is unstable against adiabatic radial per- 
turbations if and only if its radius R and geometrized mass (2M - Schwarzschild 
radius) satisfy 

r 1 - 4/3 - K(2M/R) • (10) 

Here K is a constant that depends on the star's structure and upon the PPN 
parameters. If K is positive, there is a relativistic instability. If K 
is negative, relativity stabilizes the star. In the Newtonian limit K = 0. 

In GR and BDT K is positive and of order unity. Values of K for poly tropic 
gas spheres, as evaluated in §IV of this paper, are tabulated in Tables 1 
and 2. 

Table 1 lists values of K for polytropic stars in the case of conser- 
vative theories of gravity. From the positive signs of the minimum values 
of K (column 3), we have the following conclusion: for conservative theories 
which are compatible with current solar-system experiments, relativistic 
corrections to Newtonian theory will always induce dynamical instabilities. 

It is interesting to note that 7 has a positive contribution to K while f3 
has a negative contribution; the same is true for the perihelion advance. 

This, together with the positivity of K, confirms the heuristic argument 
given at the beginning of this section. 

Table 2 lists the values of K for the general PPN formalism and for 
several particular non- conservative theories. The third column gives mini- 
mum values of K corresponding to current experimental limits on the PPN 
parameters. If £ or (which are undetermined by experiments to date) 
were sufficiently negative, then K would be negative. For example, for 
the currently viable cases = 0.76, P = 1.34, = - 0.5, C 3 = ^4 = ’ 

and {7 = P = 1, £ 2 = - 2.2, £ 3 = ^ = oj the value of K is negative. 
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Therefore we arrive at the following conclusion due to the lack of experi- 
mental information on and f, , it is inconclusive whether relativistic 

^ d 4 

effects will actually stabilize or destabilize stars. From the last three 
columns, one may notice that the Vector-Metric theory (Will and Nordtvedt 
1972) and the Papapetrou (1954a, b, c) theories have the same K-values as 
general relativity, while K-values for the Modified Yilmaz theory (Ni 1972) 
are all negative. 

Other aspects of dynamical stellar pulsations are also investigated 
in this paper. The Schwarzschild criterion is found to hold for the onset 
of dynamical instability against non-radial oscillations (convection). 
Sufficient conditions for self-adjointness of the linearized pulsation 
equations are derived. These conditions together with the condition ^ = 0 
coincide with Will's conditions for the existence of ten post- Newtonian 
conserved integrals. 

In §11 the PPN formalism is summarized, the linearized pulsation 
equations are derived, and "preferred- frame terms" (which lead to vibrational- 
secular and other Machian-type instabilities) are separated out of the pul- 
sation equations and reserved for study in a future paper. Section III 
derives a variational principle for dynamical stellar stability. Section 
IV derives the post-Newtonian conditions for the onset of a dynamical insta- 
bility. Section V derives the Schwarzschild criterion for non-radial insta- 
bilities. Concluding remarks are make in §VI. An Appendix treats the 
problem of self-adjointness. 

Throughout this paper, we follow closely the methods of Chandrasekhar 
(1965b), and we use geometrized units. The notations and conventions of 
this paper are the same as those of Chandrasekhar (1965b), and Will and 
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Nordtvedt (1972) — unless otherwise specified. 


II. PPN FORMALISM AND EQUATIONS OF MOTION 
FOR SMALL OSCULLATIONS ABOUT EQUILIBRIUM 

In the PPN formalism one describes the response of matter to gravity 
by the "local law of energy-momentum conservation" 

v . T = o (11) 

(where X is the stress-energy tensor, and V is the covariant derivative 
with respect to the PPN metric); and one describes the generation of gravity 
by matter in terms of the PPN metric (Will and Nordtvedt 1972): 

g Q0 = 1 - 2U + 2PU 2 - (2 7 + 2 + + ^)$ 1 + ^ O 

- 2[ (37 - 2(3 + 1 + $ 2 )$ 2 + (1 + C 3 )^ 3 + 2(7 + ^)$^] 

+ (a x - a 2 - a 3 )w 2 u + a g w a w P u^ - (2a^ - a 1 )w a v a , (12) 

g 0a = 2^ 7 + 3 + a i - a 2 + ^1 ^ V q: + 2^ + a 2 " 

+ - 2a 2 )w a u + a 2 w P u^ , 

gap = “ (1 + 27U)5 ap * 

Here w is the velocity of the chosen coordinate frame relative to the "pre- 
ferred-frame" of the Universe (if any); and 

. P(x', t) 

U( ~ ,t} = -f T; ~ ? T d ~' ’ (13) 
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$l(x, t) = J 


p(x', t)v 2 (x', t) 


£ - x 


, I dx' , 


5 >o( x , *0 = J 


$,(x, t) = J* 


p(x’, t)U(x', t) 


X - X 


•I dx' , 


p(x', t) n(x', t) 


X - X 


. 1 dx ' , 


p(x', t) 


H 1 **) ° j' x'| d S' < 

p(x',t) [y(x',t) • (x - x ')] 2 
tf(x, t) = J‘ — 2 = 1 dx* , 

~ « i , iO 


X - X 


Vs-*) = J 


,, p(x', t) V (x', t) 


or 

lx - x' 


dx' , 


(lb) 

(15) 

( 16 ) 

(17) 

(18) 
(19) 


W Q (x, t) . J 


p(x',t) v (x',t) (x - X ') ( 


e "B - '*a - V> 


|x - X 


. .3 


u ag(;> c > = J - 


P (x',t) (x a - x a ')tx g - y) 


x - x 


. .3 


dx' 


dx' , (20) 


( 21 ) 


Here p is rest-mass density, p is pressure, and H is specific internal 

OL 

energy all measured in the matter's rest frame, and = dx /dt is the 
matter's coordinate velocity. 

The equations of hydrodynamics governing a perfect fluid follow from 
equations (ll), (12) and the form of the stress-energy tensor (Will 1972a): 

^ ( a v “) + -^ (««V) - P + ^j ptl + (37 - 1)U1 | 

+ P ^[(57 - 1) u v a - \ (by + ^ + o L ) V a - \ a x U w«l . (22) 


8 



- |(« 2 - + 1) P ^(W a - V Q ) 


dV. 


+ p[ (^7 + ^ + o^) V p + (a x - 2a 3 ) W P ] — | 


8x 


P a [^ ■ 2 ^1 2 a 2 w ^ U q£ + a 2 W X ’ po^ 


Bx 


- P [20 - -| w • v + (a 2 + a - a^) w 2 ] 

dx 


= 0 


and 


where 


a 0 * a(p*v a ) 

a* a 


(22 cont' d. ) 


(23) 


a=p(i+v +2U + 1I + p/p) 


(24) 


0 = l(a 3 + 27 + 2 + 5 t ) V 2 +~(3 7 + 1 + - 20) u +|(1 + 5 3 ) n +|(7+? 4 )p/p (25) 


0 = \(y + 1) v 2 + - 5(37 - 2 P + 1 ) u + I n + I 7 p/p 


P* = p(l + \ V 2 + 3/U) 


(26) 

(27) 


X(x, t) = J‘ p(x', t) |x - x' | dx' 
p(x', t) 0(x', t) 

= I — — d *’ 


(28) 

(29) 


Consider an equilibrium spherically symmetric distribution of matter. 
The equation of hydrostatic equilibrium, which follows from (24), is 


d 

dr 


[1 + (37 - 1) U] p 


p[ 1 + (3 a 2 + i a 3 - 1 V ” 2] If + 2p (If + 0 If) • (30) 


Let the equilibrium configuration be slightly perturbed and describe its 
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perturbation by a Lagrangian displacement of the form 


6(x) 


ifit 


( 31 ) 


The linearized form of the equations governing the perturbation, as derived 
by combining equations (22), (23), and (30), is 

of + (a 2 - ^ + l) p(v a - w a ) + [ (57 - 1) - \(ot x + 4 7 + 4) v a ] 


+ in 


dv r 


i a t»u» a - ip(aj - 3^) J 7-I + a 2 P ^ | a^o 

dx dx dx 


dU 


dx 


a 


dx 


a 


[1 + (37 - 1) u] Ap + (37 - 1) pAU J - — [ 1 + |(« 2 + « 3 - C^) w 2 

p dx 


+ (37 - 1) U]p 


- rfl + + l°3 -I'V " 2 > T5- + 

dx dx dx dx 


+ I pa 2” 7 ” 8 T ~5 (a \ 6 - 5 5 75 ^ ’ 


(32) 


and 


AP = - P div | 


(33) 


Here and henceforth a, 0, and p* are defined by equations (24), (25), 
(26), and (27), with all quadratic velocity terms (v terms) omitted; and 

V and W are given by definitions (19) and (20) with v (x* ) replaced by 

CC OC ~ 

l a (x’)- [New V and W equal to (l/ift) times old V and W .1 The symbol 

V ~ /L a a o. a 

A denotes the Lagrangian change in the quantity that it qualifies. 

Now we must evaluate the Lagrangian changes for various quantities 
explicitly in terms of From the definition of p* and equation (33), it 
follows that 

Ap = - p(div 5 + 37AU) , (34) 
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correct to post-Newtonian order. Similarly, the first law of thermodynamics 
and the definition of "the adiabatic index" lead to 


2 

ah = (p/p ) • ap , ap = r^p/p) ap 


(35) 


respectively. Therefore, 


A0 = |(3 7 + 1 + £ 2 - 23) AU - -|^[(3r i - l )( 7 + ^) + 1 + ^)](div | + 3 7 AU) (36) 


and 


A0 = |(3 7 - 2p + 1) AU - | £ (3r i7 - 3 7 + l)(div | + 3 7 AU) . (37) 


Finally, the expressions for AU, AU , A(, and AI> can be written down from 

75 

equations (13), (21), (28), (29) as follows: 


AU = | • VU + J* p(x') r(x') 

~ ~ V ~ dx' 

p(x' ) AU(x' ) 

- 3r 4 s’! d ?' 


a' x - x' 


dx' 


(38) 


dx 


„ N (x - X ' ) ( X - X ' ) 

AU = 5 • VU + J! p(x') 5 a (x') — | T • — 2 2 ^ — dx' 

7 0 ~ 7b ~ ~ ~ *„ a i„ „i 1 5 


x - x 


- 3 7 r p(x' ) AU(x' ) 
V 


(x - X ' ) (x -x') 

7 7 5 5 


X - X 


. |3 


dx' 


(39) 


AX = | * + J p(x') | a (x') — |x - x’ | dx' 

V dx 


and 


- 3 7 J p(x' ) AU(x' ) |x - x' ] dx' 
V 


/» = i • to + J v p(*') 0(;') i a (x') ~ x ~ x 

dx ~ ” ~ 


rr d ?’ 


J, 


p(x' ) A0(x' ) 


V X - X 


dx' . 


(*0) 


(M) 
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The last two terms in equations (38), (39), (**■()), and (hi) make up the 

Eulerian changes in U, U , X and $, corresponding to their repsective 

75 

Lagrangian changes. 

Notice that the linearized pulsation equation (32) is not invariant 
under rotations. Terms linear in w couple "i-modes" (modes with spherical- 
harmonic index l) to (£ - l) and (l + l) modes; terms quadratic in w couple 
i-modes to (f - 2), (f - l), (i + l), and (i + 2) modes; all other terms 
are invariant under rotation. The terms linear in w have imaginary coeffi- 
cients; therefore they (like viscosity, energy generation, and radiative 
transport) contribute to the vibrational-secular stability of the star, but 
do not affect its dynamical stability. Terms quadratic in w contribute to 
the dynamical stability and couple different angular modes. We will delay 
until a later paper all analyses of w-dependent terms ("preferred- frame 
terms") — including both the problem of vibrational-secular stability (linear 
in w) and preferred- frame influence on dynamical stability (quadratic in 
w). Thus, we shall set w = 0 throughout this paper. 

III. THE VARIATIONAL PRINCIPLE 

Equation (32), when supplemented by the expressions for the various 
Lagrangian changes in terms of £, becomes explicitly an equation for £. As 
boundary conditions, we shall demand that Ap = 0 at the surface of the star 
( r = R), and that there be no physical singularity at the star's center 
(r = 0). Equation (32) together with the boundary conditions then con- 
stitutes a characteristic value problem for ft. 

If a characteristic value problem is self-adjoint, the orthogonality 
relations for its characteristic functions hold and a variational base for 

(' 
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determining Q, can be obtained. The stellar pulsation equations in the post- 
Newtonian limits of general relativity and Brans-Dicke theory are self-adjoint 
(Chandrasekhar 1956; Nutku 1969) as is the equation for radial oscillations 
in the full theory of general relativity (Chandrasekhar 1964b). In the 
Appendix, it is shown that the characteristic value problem in the PPN for- 
malism is self-adjoint if and only if 


“l = a 2 = a 3 = ^2 = ^3 = k = ° • (1|2) 

Although, in the general case, the characteristic value problem is not self- 
adjoint and the orthogonality relations do not hold, a variational integral 
can still be constructed in the following manner: 

Take equation (32) with Q, replaced by the characteristic value for 

the i-th normal mode, with g replaced by the corresponding characteristic 
function g^, and with w-terms deleted. Dot into this equation g^, the 
characteristic function for the j-th mode, and integrate over the interior 
of the star. Thereby obtain 




(1+3) 


(i i) ' 

Here Q v is expression (A. 4) with 1 replaced by i, 2 replaced by j, and 

complex conjugations deleted; ^ is the symmetric part of the right-hand 

side of (A. 3), with similar replacements; and ^ is (A. 4) with similar 

replacements. Notice that R^' ^ is of post-Newtonian order: 


^ = 0 ( 2 ) . 


(44) v 


^By 0(n) we mean, in Chandrasekhar's (1965a) language, 0(c n ), 
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From equations (43) and (44), it follows immediately that the standard 
orthogonality relation for characteristic functions is valid to Newtonian 


order, i.e. 


jC 1 , j) = 0 (2) / n^] . 


(45) 


We assume, without proof, that the characteristic functions || 
complete set; and we normalize them to give 


L) 


from a 


,(I,D = ! 


(•* 6 ) 


(no summation on capital letters). 

Let be a solution which differs from g^^ by post-Newtonian order 


and has norm 1, i.e. 


^(i) = |( I ) + 0 (2) 


q( p ? (i) , p s ( i) ) ‘ 1 • 


Expand ^g^ in terms of g^^: 


P I (I) = s c xj | (j) , 

and from equations (45), (46), (47), and (48), obtain 


(47) 

(48) 


c zj = l + 0(4), (j = I) 


C TJ = 0(2), (j/1) . 


(49) 

(50) . 


By combining equations (43), (45), (49), and (50), we obtain 

s(V*>, ',(«) t P S (I) ) - .z U (J) ] 2 c tJ C tj , 

Jj J 


= [n (l ^] 2 Q (I,J ^ + 0(h) ; (51) 
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by combining equations (48), (4-9), and (50), we obtain 

Q ( P | (I) , P | (I) ) = C Ijt Q (j#j,) = Q (l,l) + 0(4) ; (52) 

■)> J 

and by combining equations (51), (42), we finally conclude that 

U (I >] 2 q( P { (I) , P s (1) ) - V 1 ’) + V 1 ') + 0(4) . 

P (i) 

Therefore we can use this equation, and any functions | that agree with 
only at Newtonian order, to calculate [n^^] 2 to post-Newtonian order. 
The Newtonian proper solutions 


( 53 ) 


inserting from 


are one set of such functions. 

I ~ 1 

By suppressing the prefix "p" and superscript labels, by 
Appendix the values of Q, S, and R, and by performing some reductions, we 
bring our variational expression (53) into the form 


Qn 2 = J* T. p[l + (3r - l) U ] ( div |) 2 dx + J* (div |) | 

TT "** ~ ~ TT ~ 


a d 


v 


dx 


a 


[1 + (37 - l)U]p 


dx 


+ J* ( 3 r ,7 -37 + 1) PAU div £ dx + J ( 2|3 - l)AU£ a -^- dx 
V ~ ~ V dx 


+ J* (3r,7 - 37 + 1) P div £ ^ dx - J* p £ 


.a cAU 


_ r a ol cAU . _ r* .u ca® , 

- 2 J p 0 | — dx - 2 J p| -SSL dx . 
V dx V dx 


dx 

a 


v 


dx 


a 


dx 


(54) 


Here 


Q = J* CT|S | 2 dx + \(a. - ^ + 1) J* J p(x) p(x' ) - 

V TT ~ ~ CL C. L TT T7 ~ 


!(x) • l(x') 


V V 


x - x 


dx dx' 


- ^i + l) L L p W p( ?') 


U(x) • (x - x')][£(x') • (x - x')] 


V V 


X - X 


. |3 


dx dx' 


«(*) 


+ (57 - 1) J* J* p(x) p(x' ) - , dx dx' - -x(4 7 + 4 + OC ) J J 

J v v 15 " 5 ~ ~ ^ vv 


p(x) p(x' ) 


l(x) • S(x') 


dx dx' 


(55) 
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is a positive definite quantity in the post-Newtonian limit since the domi- 

nant, Newtonian part, J v a||| dx, is positive definite. This equation can 

2 

be used for a variational determination of Q . 

We shall now analyze the Lagrangian displacement £ into normal modes 
belonging to different vector spherical harmonics. (Since the pulsation 
equation without w-dependent terms is invariant under rotation, this pro- 
cedure is justified.) Following the procedure of Chandrasekhar (1961, 1964c) 
and Lebovitz (1965), we define: 


and 



(56) 

(57) 







(58) 


(5 , £„, 5m are physical components, not covariant components). After 
r 0 F 

manipulations similar to those in §IV of Chandrasekhar (1965b), we obtain 

2 

the following expression for the variational determination of ft : 

v? - j E r x p[i + (37 - i)u]|-£ (+ - *)J 2 ^§ 

+ 2 X* Tr u + & - 

- Wrr I'd + 2 * ♦ s*> W, -5F - K, *?)■*' 

+ (28 - 1 - £ ? ) J* R p[AU(r)] 2 r 2 dr 
0 

+ (6r i 7 - 67 + 2 + 3r x ^ - 3£ 4 + ^) J* P AU(r) (* - X) dr . (59) 
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where 


v*>-r p<°> 4 ^**£4* - 


K. 


(r) = J R £M U + D^-^ 


r s 


s ds 


ds 


and 


AU(r) = 


bn l J l< r ’ 


bl 


«_|i 

+ 1 1 

*• r 


. , - r l K (r) 

f+1 ' 


^ dU 

+ “2 cf r 
r 


( 60 ) 

( 61 ) 

( 62 ) 


IV. THE POST- NEWTONIAN CONDITION FOR THE ONSET 
OF DYNAMICAL INSTABILITY 


Consider the case of radial pulsations, i.e. pulsations with 

l = 0 and X = 0 . (63) 

The substitutions 

^ = r 3 ^ and l r = r T] , (6h) 

reduce equation (59) to the form 

Qfi 2 = J 1 p[ 1 + (37 - 1) U ]| Fl T h (^}) 2 + (3I-J - !•) ~ (r% 2 )]dr 

+ (2P - 1 - 5 ) J R p[AU(r)] 2 r 2 dr 
0 

+ (6Tj7 - 67 + 2 + 3 r^ u 

- < 3r A ' 3 ^ + ^ J 0 R 
+ ^2 J* Q AUr l r dr ’ 

where 


- + c 3 ) J* pAU (r ) ( r tj) dr 

dU d(qr 3 ) 
p d7 ^ dr r dr 


( 65 ) 
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( 66 ) 


AU(r) = - bn J* R p(s) ^(s) \ ^ 

r - 


s r 

Recall that p and p are the distributions of pressure and density in the 
equilibrium configuration in the post-Newtonian approximation, and they 
therefore include terms of 0(2). 

The condition for marginal stability (instability) follows from equa- 

2 

tion (65) by setting (1=0. In the particular case P^ = const. — which 
implies q = const, at the point of onset of instability in the Newtonian 
limit, i.e. q = const. + 0(2) — the condition for marginal instability (eq. 
[65] with ft = 0) involves the structure of the equilibrium configuration in 
the Newtonian approximation alone. Under these conditions the criterion 
for marginal instability becomes P^ = T cr ^ t , where 


crit 


Here 


k JL_ 
3 + 3W 


(2P - 1 - £ 2 ) J* R [AU(r)] 2 dM(r) + 3[2y + 2 + ^ 


J R f AU(r) dM(r) 


- 3 <?3 + /* l (f ■ ') + t 8 S f * (|?) • (57) 


W = - 12* J* R pr 2 dr , 
0 

. P 

dM = 4npr dr , 


( 68 ) 

(69) 


and 


AU(r) = - 4rt J R psds + r 


dr ' 


(70) 


This result agrees with those in general relativity (Chandrasekhar 1965b) 
and in Brans-Dicke- Jordan theory (Nutku 1969) when specialized to the cor- 
responding PPN values. 
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Criterion (67) for marginal instability may be reformulated as follows. 
A dynamical instability will set in if and only if the following inequality 
is satisfied: 


_ _ k __ 2M 

r crit 3 + K R 


(71) 


Here M is the mass and R is the radius of the configuration, and K is a 
constant (typically of order unity), depending on the Newtonian structure 
of the configuration. If K is positive, there is a relativistic instability; 
if K is negative, then relativistic effects stabilize the star. 

For polytropes, an explicit expression for K can be obtained, from 
equations (67)-(7l), in terms of Lane-Emden functions: 


K = | (1 + 7 ) \ - (2p - 1) K g + (5 3 + £ u )({ K x - 3K^) + ^(^ - K^) (72) 


where 


*1 = 


= 


*3 


- n) |» llg + . iq i ||e n +l *2 ^ 

X6(n + 1) t V,'! 3 ° I 11 J ? 1 


r~773 ' £ [ e * hlVl] 2 ®” “6 


(5 - n) 
181 


1 1 1 

5 - n 


18(n + 1) 1^10^ j 3 '0 


J 1 ® n+1 dft 


(73) 


W 


(75) 


\ ■ 18 y ,3 • ft® + ,‘iivi] ^ |3 <76) 


and where n is the poly tropic index, ^ is the first zero of the Lane-Emden 

function 0 , and ©, ' is the value of the derivative of © at L. The 
n' 1 n 1 

values of the constant K, evaluated with the aid of the foregoing formula 
for various values of n, are listed in Table 1 for conservative theories 
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and in Table 2 for the general PPN formalism and for non- conservative theories 
See §1 for discussions of these tables and for the conclusions inferred from 
these tables. 


V. NON-RADIAL OSCILLATIONS AND THE SCHWARZSCHILD 
CRITERION FOR CONVECTION 


We shall now obtain the condition for the occurrence of a neutral mode 
of non-radial oscillation belonging to I > 1 in the general PPN framework. 

As in the last two sections, all "preferred- frame" effects (w-dependent 

2 

terms) will be ignored. By setting ft = 0 in equation (32), by following 
an analysis parallel to §VI of Chandrasekhar (1965b), and by using the 
result of §VIII of Chandrasekhar (1965b), one can show that, to 0(2), a 
necessary and sufficient condition for the occurrence of a neutral mode of 
non-radial oscillation is that 


S(r ) 1 + n 


P 3 ' P 


k p dp/dr 
- 1 p dp /dr 


= 0 


(77) 


i.e. that 


S(r) = 0 


(78) 


over some finite interval of r. Here 


S(r) 


i*p _ r p i£ 

dr 1 1 p dr 


(79) 


is the "Schwarzschild discriminant", and a nd are defined by 


and 


1 + 


d(log T) 
d(log p ) J s 




n = 



p 

p 


(80) 


(81) 
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By following a procedure similar to §VII of Chandrasekhar (1965b), one can 
also derive this condition from the variational principle (61 ). 

The proportionality of the Newtonian and the post-Newtonian discrimi- 
nants implies that the physical condition for convective instability remains 
the same in the PPN formalism as in general relativity and in Newtonian 
theory. Although for some PPN values, the characteristic value problem is 
not self-adjoint, the Schwarzschild criterion is still valid, and no new 
dynamical instabilities occur. 

VI. CONCLUSIONS 

In this paper, stability criteria for stellar pulsations were found 
using the general PPN formalism. These criteria are valid for almost all 
metric theories of gravity in the post- Newtonian opproximation, when one 
ignores preferred- frame effects. As in general relativity, so also for 
conservative theories (conservative theories do not have preferred- frame 
effects), the relativistic corrections do actually induce dynamical insta- 
bilities in stars. But in the general case the present experimental uncer- 
tainty in the PPN parameters | , makes it inconclusive 

whether relativistic effects will actually stabilize or destabilize stars. 

As experimental tests are performed to higher precision, the answer may 
become definite. The differences in the dynamical stability criterion for 
various theories may affect the evolution of white dwarfs and supermassive 
stars; such effects are worth exploration. The relationship of non-self- 
adjointness to the non-existence of conservation laws is intriguing and 
should be examined further. 

A subsequent paper will deal with the problem of Machian instabilities 
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due to preferred- frame effects (w-dependent terms). Those instabilities, 
when combined with astronomical observations on white-dwarf pulsations, may 
lead to tight experimental limits on the "preferred- frame parameters" 0!^, 
Qv,, an d 0^ * 
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APPENDIX 


SELF-ADJOINTNESS OF THE CHARACTERISTIC VALUE PROBLEM 


We shall here derive the constraints which the PPN parameters must 
satisfy for the characteristic value problem [ eq. (32)] to be self-adjoint. 
For this prupose we do not delete the w-dependent terms ab initio (cf. end 
of §11). 


By bringing the right-hand side over to the left, write equation (32) 
in the form = 0, where ^ is a linear operator. This equation is self- 
adjoint [or can be made so by multiplication with some weighting function 
y(x)] if and only if 


J* I 


(D* . 


y =£ i 


( 2 ) 


d ? = J* b 


;[ 


( 2)" 1 


* y*£ 6 


(1) 


dx 


(a. l) 


where the complex conjugation is not permitted to act on the eigenvalue 
Q (which is contained in=£). In this equation and are arbitrary 

functions satisfying the boundary conditions at r = 0 and r = R (but not 

necessarily satisfying £ | ^ ^ = 0 or a? = 0); V is the interior of the 

12 3 

star; denotes complex conjugation; and dx denotes dx dy dz = dx dx dx • 

From this definition one readily verifies that (i) if the weighting function 

y is chosen real, then the in- terms prevent self-adjointness; (ii) if y has 

2 

any imaginary part, then the n -terms prevent self-adjointness. It is pos- 
sible to get rid of the in- terms by demanding = 0; but it is 

2 

not possible to get rid of the n -terms. Therefore, to have any hope of 
self-adjointness one must choose y.(x) real and 


0^ = a 2 = 0^ = 0 . (A. 2) 

Insist, then, that = a = 0; and try, for the moment to prove 
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self-adjointness with the trivial weighting function y(x) = 1. Then the 
terms on the left-hand side of equation (32) give, when integrated, 


[n]‘ 


J* cr |^ 2 ) a d x + J* J p(x) p(x' ) 

V V V 


e (l)a*, v e (2 )oc, ,s 
I (x) | v (x') 


X - X' 


dx dx' 


- ¥ a 2 - ?i + l) X J 1 oM p(;') ! (1) "(;)(*“ - *'“) (tetyz'HiP - 


X 


,3 


V V 


dx dx' 


X - X 


.,3 


(5r - 1 ) J J* p(x) p(x') ^ 1 ^ a *(x) |^ 2 ^ a (x) 


V V 


dx dx' 


x - x 


dx dx 1 I 

^(a x + ky + k) p(x) p(x') ! (l)a *(x) ^ 2 ^ a (x') =n 2 Q (1,g 


[For use in the text of the paper we retain all a terms except those that 

depend on w; but we keep in mind that henceforth in this appendix the a' s 

(l 2) 

vanish. ] Notice that Q ’ is manifestly self-adjoint in the superscripts 

(l 2) (2 l)* 

1 and 2, i.e. Q v ’ ' = Q v ’ . Reduction of the right-hand side of the 

integrated equation is less straightforward; it requires integrations by 
parts, followed by substitutions for various Lagrangian changes, simplifica- 
tions, and rearrangements. After some work, the right-hand side is brought 
into the form 


fi 2q(!,2) = (• r n + (37 - 1)U] p div £^* div dx 

V TT 1 ~ ~ ~ 


r A jL 

v„ a dr dr 


V 


[1 + (37 - 1 ) U ] p 


[x • | (l) *][x • | (2) ] 


dx 


(A.4) 


2k 


. (A.3) 



dx 


rx • x • p (2) 1 

+ [1 + (37 - 1)U] |£ div e (2) +- f div 

~ f I I 1 + 2[?(S) + ?(*')]} 

V V 1 1 


»(l) 


.( 2 ) 


[Ap V (x) - (x) • Vp(x)][Ap (x’) -£ 

" lx - x' I 


+ J*„ P Ji + 20) 


x • ? (1) * x • f (2) 1 

2 « ( 2 ) * 5 (i)J 

AIT ' + AIT ' dx 

r r 1 ~ 




(!)* x . .(2) / 

(2) ~ I .Jl)*l 


An' 7 + 


ATI |dx 


- (37r, -37 + 1) f P[AU^^* AH^ + Al/ 2 ^ AH (l) *]dx 

l v 


- 2 


p [dU dp 
J v |dr dr 


(0 - 0 ) + 


dU d0 
dr dr 


[x • £ (l) *][x • | (2) ] 


dx 


- (67 + 1 + - 2 P) J PAU^ 1 ^* AU^ 2 ^ dx + R^ 1 ’ 2 ^ 

d V 


(A. k 


where 


t (l,2) = [- 3^(1^ - 1) - 5 3 ] J* pAIT (l) * AU (2) dx 


( 2 ) 

. .( 2 ), 


+ 2 J / d fW ‘« 


(2) dU 
+ ^ d? 


« ; 

- ait 

AU^ 1 ^* Ap^ 2 ^ 

+ AU 

[x • | (l) ] ] 

dx 

r J 


.( 2 ). 


dr 


(x‘) * Vp(x')] 

— ~ — dxdx ' 


cont' d. ) 


(A. 5) 
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Equation (A.k) reduces to equation (A. 6) of Chandrasekhar (1965b), if we 
substitute in the PPN parameter values for general relativity. Aside from 
the terms on the right-hand side of equation (A.k), like of 

the left-hand side, are manifestly self-adjoint in the superscripts 1 and 2 . 
Therefore, the condition for equation (32) to be self-adjoint with weighting 
function 5 /.(x) = 1 is 

r (1,2) = R ( 2, 1 )* (A. 6 ) 


or, equivalently (since one insists on this equality for all choices of 

? 2 = = ?!, = 0 ■ ( A - 7 > 


and |^ 2 k: 


Might some other choice of weighting function aside from constant per- 
mit one to relax these constraints and still retain self-adjointness? No; 

o 

because any other choice of y(x) will destroy the self-adjointness of the 

2 (1 p) p 

left-hand side (fi Q v ’ ), and the arbitrariness of fi will prevent the 

non- self-adjoint terms thus created from always cancelling non-self-adjoint 

terms on the right-hand side. 

^3 = 

= 0 — are necessary and sufficient conditions for the self-adjointness 
of the linearized pulsation problem in the PPN formalism. These conditions, 
together with the condition that ^ = 0, are just Will's (1971b) conditions 
for theories of gravity to have post-Newtonian conserved integrals for 
energy, momentum, angular momentum, and center- of-mass motion. 


In summary, equations (A. 2) and (A. 7) — i.e. 0!^ = Ctg = <X^ = 


C 
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angular momentum, and center-of-mass motion — i.e. a theory with Ofy = ag = QSj = £i = = £3 = = ()• 

These are minimum values compatible with current experimental limits on (3 and 7: p < 1.34, 7 > 0.76. 

The values in parentheses are those derived by Chandrasekhar (1965). 

The values in parentheses are those derived by Nutku (1969); the discrepancies are presumably due to errors 
in Nutku' s numerical calculations. 











VALUES OF THE CONSTANT K FOR POLYTROPIC STARS 



28 



FOOTNOTES FOR TABLE 2 


These are minimum values compatible with current experimental limits on (3, 

7, and £ ; p < 1.3^, 7 > 0.76, > - 0.05. 

^Will and Nordtvedt (1972). 

c Papapetrou (l95^a, b, c). 

^Ni (1972). The values of K depend on which "matter density" one chooses as source 
for the gravitational field: p = T^ u* u^ = component of stress-energy tensor 
along four-velocity of matter [upper values]; or p = trace ( T^, ) [lower values], 
(cf. Ni 1972.) 
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